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Abstract. This note illustrates how a simple random walk with possibly long jumps is related to fractional 
powers of the Laplace operator. 
', The exposition is elementary and self-contained. 
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' The purpose of this note, which is mainly pedagogical, is to show in a simple, concrete example how 

p ^ I singular integrals naturally arise as a continuous limit of discrete, long jump random walks, and to recall 

a simple description of the integral kernels in terms of the Fourier multipliers. 

Singular integrals and nonlocal (especially fractional) operators are a classical topic in harmonic analysis 
^ • and operator theory [Lan72[ ISte70| and they are now becoming impressively fashionable because of their 

connection with many real-world phenomena. 

Indeed, nonlocal operators arise in the thin obstacle problem |Caf79) . in optimization |DL76| . in fi- 
nance |CT04| . in ph ase tra nsitions |AB98l lABSOSl lCSSTOSl ISVOSbj . in stra tified m aterials |SV08a| . in 
anomalous diffusion [MKOO] . in crystal dislocation [Tol97j . in soft thin films [Kur06| . in some models of 
semipermeable membranes and flame propagation |CRS08j . in conservation laws [BKWOl] , in the ultra- 
i relativistic limit of quantum mechanics [FdlL86], in quasi-geostrophic flows |MT96[ ICorQSj . in multiple 

scattering [DG75, CK98, GKOlj, in minimal surfaces [CRS07], in materials science Bat06 and in water 
waves pto57. .Zak68. .Whi74, ,CHS92l ICCTMl INSMI [CW95. .dlfcpMl l(]SS971 ICNHnl IGGOa/.HNnHl iNTOSl 
ldlLV08| . See also [Sil05| for further motivation. 

From a probabilistic point of view, such nonlocal operators are related to Levy processes jIto84| IBer96| 
I IBG991 IJMW05] . A naive example in Section [T] will show a possible probabilistic interpretation. Then, in 

Section [2l we will recall an easy recipe for representing the singular integral kernel in terms of the Fourier 
symbols. The particular (and particularly interesting) case of the fractional Laplacian will be discussed in 
^ . Sectional 

Jh I The exposition is self-contained and no prerequisite is needed: a basic undergraduate math knowledge 

will suffice. 

In the rest of this note, we will discard all the multiplicative normalizing constants. That is, following 
a convention typical of the lectures on Fourier analysis, we will write X — Y to mean that there is some 
normalizing constant C > sue ]E that X = CY. 

1. Long jump random walks and singular integral kernels 
Let 3C : R" ^ [0, +oo) be even, that is X{y) = OC{-y) for any y eW\ and such that 

(1) J2 ^(^) = 1- 

Give a small h > 0, we consider a random walk on the lattice ft,Z". 



It is a pleasure to acknowledge the interesting conversations on the subject of this note with Domenico Marinucci. The 
author has been supported by MIUR Project Variational Methods and Nonlinear Differential Equations and FIRB Project 
Analysis and Beyond. 

"'^We hope that no reader is bothered by the fact that this convention implies, for instance, that 2it = 1. It seems not to be 
a joke that on February 5, 1897, the House of Representatives of the State of Indiana unanimously passed a bill which would 
have supported such a new mathematical truth |Ind97) . 
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We suppose that at any unit of time t (which may depend on h), a particle jumps from any point of /iZ" 
to any other point. 

The probabihty for which a particle jumps from the point hk G /iZ" to the point hk is taken to 
be DC(fc-fc) ^X{k-k). 

Note that, differently from the standard random walk, in this process the particle may experience 
arbitrarily long jumps, though with a small probability. 

We call u{x, t) the probability that our particle lies a,t x G hl/^ at time t G rZ. 

Of course, u{x^ t + t) equals the sum of all the probabilities of the possible positions x + hk at time t 
weighted by the probability of jumping from x -\- hk to x. 
That is, 

w(x,t + r) = ^ %{k)u{x + hk,t). 
fceZ" 

Therefore, recalling the normalization in ([T]), 

(2) u{x, t + r)- u{x, t) ^ %{k) {u{x + hk, t) — u{x, i)) . 

Particularly nice asymptotics are obtained in the case in which t — h°' and 3C is a homogeneous kernel, 
say, up to normalization factors 

(3) D<:(y) = |yr (for y ^ and, say, 3C(0) = 0), 
with 

a G (0,2). 

We observe that ((T]) holds (again, up to normalization) and 

(4) = h'^Xihk). 

T 

Thus, in this case it is convenient to define 

i'iy, X, t) = X{y) {u{x + y,t) - u{x, t)) 

and to use (0]) to write as 

u{x,t + T) -u{x,t) ^ 3C{k), ( , ui +\ ( 
= [u(x + hk,t) — u[x,tjj 

(5) = /i" ^ 3C(/ifc)(u(a; + hk, i) — u(x, <)) 

= /i" ^ i){hk,x,t). 

Since the latter is just the approximating Riemann sum of 

i:{y,x,t) dy, 

by sending r = /i" — > 0^ in ([51), that is, by taking the continuous limit of the discrete random walk, we 
obtain 

dtu{x,t) / ^{y,x,t)dy 

that is 

ir\ p, ( +\ f u{x + y,t)-u{x,t) 

This shows that a simple random walk with possibly long jumps produces, in the limit, a singular integral 
with a homogeneous kernel. 
We remark that the integral 

[ u{x + y) -u{x) 
which appears in ^ has a singularity when y = 0. 
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However, when a G (0, 2) and u is smooth and bounded, such integral is well defined as principal value, 
that is as 

hm / "(^ + y) -"(^) 
Indeed, \y\~^'^^"^ is integrable at infinity and 

Vulx) ■ y 



= 



as principal value, because the function y/|2/|"^" is odd. 

Therefore, we may write the singular integral in ^ as principal value near in the form 

u{x + y) — u{x) — 'Vu{x) ■ y 



L 



■ dy 



and the latter is a convergent integral near because 

\u{x + y) - u{x) - Vu{x) ■ y\ ^ \\D'^u\\l'^ |yp 



|y|n+Q |y|ri+Q! |y|ri — 2+Q 

which is integrable near 0. 

It is also interesting to write the singular integral in ([7]) as a weighted second order differential quotient. 
For this, we observe that, substituting y = —y, we have that the integral in ^ equals to 

(8) / "(^ " dy. 

Therefore, relabeling y as y in ([8]) , we have that 

u{x + y) ~ u{x) 

V\ 

u{x + y) ~ u{x) f u{x — y) — u{x) 

u(x + y) + u{x — y) — 2u{x) 



2 / ^ — -4^- — ^dy 

(9) = / ; ^^ ""^-^ dy+ I "^"^ , "^"^^ dy 



■ dy. 



The equality obtained in ^ shows that the singular integral in ([7]) may be written, up to a factor 2, as an 
average of the second incremental quotient u{x + y) + u{x — y) — 2u{x) against the weight 

Such a representation is also useful to remove the singularity of the integral at 0, since, for smooth u, a 
second order Taylor expansion gives that 

u{x + y) + u{x - y) - 2u{x) ^ \\D^u\\l'^ 



^ |y|n— 2+Q 

which is integrable near 0. 

It is known [Sil05" that the singular integral in ([7]) is related to the fractional Laplacian (— A)"/^. This 
relation will be outlined here below (see, in particular and below). 

It is also interesting to write the displacement of the above random walk at time nr, for any n G N. 
Namely, if hej € /iZ" is the jump performed at time jr (that is, the "innovation"), the above discussed 
random walk is made in such a way that the probability that £j equals k is DC(fc). 

The displacement at time nr is then the sum of these innovations, that is 

n 

The /3-moment associated to this process is then 

fceZ" fcGZ"\{o} 

which is finite if and only if /3 < a. 

In the probability theory framework, this is interpreted as the innovation being in the domain of attrac- 
tion of an "a-stable random variable" |ST94| . 

In particular, the associated variance is not finite, thus refiecting that the process is not Gaussian. 
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2. Kernels and Fourier symbols 

Given a "nice" (say, smooth and with fast decay, for simphcity) function u, the long jump random walk 
of Section [1] has lead us to the study of integrals of the type 

(10) / {uix + y)+uix~y)^2uix))Xiy)dy, 
due to dHl). 

If we call £jU the integral in (fTO|) . one may consider L a linear operator and look for its "symbol" (or 
"multiplier") in Fourier space. 

That is, if 3^ denotes the Fourier transform, one may think to write 

(11) Lu{x) ^3'-^{§ {S'u)), 

for some function § : R" M. 

The interesting fact is that 3C and § are related as follows: 

(12) m= I (cos(e-2/)-l)3<:(y)dy, 

up to normalization factors. 

To check that ([12]) holds, one simply Fourier transforms pT|) in the variable calling ^ the corresponding 
frequency variable: making use of pO|) one obtains 



§(0(^")(e) = ^(M 

= {u{x + y) + u{x-y)-2u{x))X{y)dy^ 

{j{u{x + y)+ u{x -y)- 2u{x))^%{y) dy 

yi-y + e-^i-y -2) {Jumny)dy 
(e^i-y + e-^^-y -2) %{y)dy (Ju)(0 
(cos(e-y)-l) X{y)dy (?w)(0, 

proving (fT^ . 

3. The fractional Laplacian 

The fractional Laplacian may be naturally introduced in the Fourier space. 
Indeed, one has that 

and therefore 

-Am = J-i(|^|2(?w)). 
Thus, it is natural to define, for a G (0, 2), 

(13) (-A)"/2u^y-i(|^|"(j-^)). 

It is known |Lan72[ ISteTOj that such a fractional Laplacian may be also represented as the principal value 
of singular integral, namely 

up to normalizing constants - again, the above integral is intended in the principal value sense. 
Notice that, by one can also write as 



/2 / u{x + y) + u{x - y) - 2u{x) 



y]n+a 



dy, 



up to normalizing factors. 

We give here a simple proof of the equivalence between the definitions in and in 
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For this, we observe that, in the notation of (fTO| and (fTT|) . we may write (fT3|) and ([Till as 

m^l^r and 3C(y)--|yr 
Therefore, by (|12p . such equivalence boils down to prove that 

(15) / '-i^dy. 

To prove first observe that, if C = (Ci, ■ ■ • , Cn) ^ -l^"! '^^ have 

1-cosCi ^ ICiP ^ 1 



near C = 0, therefore 

(16) / —TT7—r^d( is finite and positive. 



We now consider the function 

We have that d is rotationally invariant, that is 

(17) a(?) = a(ieiei). 

Indeed, if n = 1, then one easily checks that d{~£,) = 3(0 j proving (fTTl) in this case. 
When n ^ 2, we consider a rotation R for which 

^(leiei) =C 

and we denote by R"'" its transpose. We obtain, via the substitution y — R^y, 

^^^^ ^ L — — 

l-cos((|g|ei).(fl^y)) ^ 

|y|n+a y 

l-cos((|C|ei)-y) ^ 

|y|„+a J' 

which proves (fT7|) . 

As a consequence of (fT7)) and the substitution C = l^ly gives that 



1 -cos(lClyi) 



|y|n+a 



1^1" A" |c/ici|"^" ' 

up to normalization factors, hence is proved, thus so is the equivalence between (|13p and (HH). 

We remark that, from (jl4p . the probability density of the limit long jump random walk in ^ may be 
written as 

(18) dtu^ -{-A)"/^u. 
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